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About 25 per cent of hot Jupiters (extrasolar Jovian-mass planets with close-in orbits) are 
actually orbiting counter to the spin direction of the star 1 . Perturbations from a distant 
binary star companion 2 3 can produce high inclinations, but cannot explain orbits that are 
retrograde with respect to the total angular momentum of the system. Such orbits in a stellar 
context can be produced through secular (that is, long term) perturbations in hierarchical 
triple-star systems. Here we report a similar application to planetary bodies, including both 
the key octupole-order effects and tidal friction, and find that it can produce hot Jupiters in 
orbits that are retrograde with respect to the total angular momentum. With distant stellar 
mass perturbers such an outcome is not possible 2 3 . With planetary perturbers the inner 
orbit's angular momentum component parallel to the total angular momentum need not be 
constant 4 . In fact, as we show here, it can even change sign, leading to a retrograde orbit. A 
brief excursion to very high eccentricity during the chaotic evolution of the inner orbit can 
then lead to rapid capture, forming a retrograde hot Jupiter. 

Despite many attempts^EHD] there is no model that can account for all the properties of the 
known hot Jupiter (HJ) systems. One model suggests that HJs formed far away from the star and 
slowly spiraled in, losing angular momentum and orbital energy to the protoplanetary disk ! 12 * 13 !. 
This "migration" process should produce planets with low orbital inclinations and eccentricities. 
However, many HJs are observed to be on orbits with high eccentricities, and misaligned with 
the spin direction of the star (as measured through the Rossiter-McLaughlin effect!^) and some 
of these (8 out of 32) even appear to be orbiting counter to the spin of the star. In a second 
model, secular perturbations from a distant binary star companion can produce increases in the 
eccentricity and inclination of a planetary orbiP^. During the evolution to high eccentricity, tidal 
dissipation near pericenter can force the planet's orbit to decay, potentially forming a misaligned 
HpE. Recently, secular chaos involving several planets has also been proposed as a way to form 
HJs on eccentric and misaligned orbits-^. A different class of models to produce a tilted orbit is 
via planet-planet scattering 5 , possibly combined with other perturbers and tidal friction 7 . In such 
models the initial configuration is a densly-packed system of planets and the final tilted orbit is a 
result of dynamical scattering among the planets, in contrast to the secular interactions we study 
here. 



In our general treatment of secular interactions between two orbiting bodies we allow for 
the magnitude and orientation of both orbital angular momenta to change (see Figure 1). The 
outer body (here either a planet or a brown-dwarf) gravitationally perturbs the inner planet on 
time scales long compared to the orbital period (i.e., we consider the secular evolution of the 
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system). We define the orientation of the inner orbit with respect to the invariable plane of the 
system (perpendicular to the total angular momentum): a prograde (retrograde) orbit has i\ < 90° 
{ii > 90°), where i\ is the inclination of the inner orbit with respect to the total angular momentum 
vector. Note that the word "retrograde" is also used in the literature to indicate orbital motion 
counter to the stellar spin. The directly observed parameter is actually the projected angle between 
the spin axis of the star and the orbital angular momentum of a HJ. Our proposed mechanism 
can produce HJs that are "retrograde" both with respect to the stellar spin and with respect to 
the total angular momentum. By contrast, a distant stellar companion can only succeed in the 
former. See the online Supplementary Information for more details; henceforth we will use the 
term "retrograde" only to indicate an orbit with %\ > 90° as define above. 



We assume a hierarchical configuration, with the outer perturber on a much wider orbit than 
the inner one. In the secular approximation the orbits may change shape and orientation but the 
semi-major axes are strictly conserved in the absence of tidal dissipatiorP^. In particular, the 
Kozai-Lidov mechanism 17 19 produces large- amplitude oscillations of the eccentricity and inclina- 
tion when the initial relative inclination between the inner and outer orbits is sufficiently large ( 
40° < z < 140°). 



We have derived the secular evolution equations to octupole order using Hamiltonian per- 
turbation theory^ESnD in contrast to previous derivations of "Kozai-type" evolution, our treatment 
allows for changes in the z-components of the orbital angular momenta (i.e., the components along 
the total angular momentum) L z l and L z 2 (see Supplementary Information). The octupole-order 
equations allow us to calculate the evolution of systems with more closely coupled orbits and with 
planetary-mass perturbers. The octupole-level terms can give rise to fluctuations in the eccen- 
tricity maxima to arbitrarily high values^H in contrast to the regular evolution in the quadrupole 
potential^SE], where the amplitude of eccentricity oscillations is constant. 



Many previous studies of secular perturbations in hierarchical triples considered a stellar- 
mass perturber, for which L Zj i is very nearly constanP^KS. Moreover, the assumption that L Zj i 
is constant has been built into previous derivation ^ 22 * 24 !. However, this assumption is only valid 
as long as L 2 3> L\, which is not the case in comparable-mass systems (e.g., with two planets). 
Unfortunately, an immediate consequence of this assumption is that an orbit that is prograde rel- 
ative to the total angular momentum always remains prograde. Figure Q] shows the evolution of a 
representative system (here without tidal effects for simplicity): the inner planet oscillates between 
prograde and retrograde orbits (with respect to the total angular momentum) as angular momentum 
flows back and forth between the two orbits. 



Previous calculations of planet migration through "Kozai cycles with tidal friction'^^^M] 
produced a slow, gradual spiral-in of the inner planet. Instead, our treatment shows that the eccen- 
tricity can occasionally reach a much higher value than in the regular "Kozai cycles" calculated 
to quadrupole order. Thus, the pericenter distance will occasionally shrink on a short time scale 
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(compared to the Kozai period), and the planet can then suddenly be tidally captured by the star. 
We propose to call this "Kozai capture." 

Kozai capture provides a new way to form HJs. If the capture happens after the inner orbit 
has flipped the HJ will appear in a retrograde orbit. This is illustrated in Figure 2. During the 
evolution of the system the inner orbit shrinks in steps (Fig. 2c) whenever the dissipation becomes 
significant, i.e., near unusually high eccentricity maxima. The inner orbit can then eventually 
become tidally circularized. This happens near the end of the evolution, on a very short time scale 
(see Fig. 2, right panels). In this final step, the inner orbit completely and quickly decouples from 
the outer perturber, and the orbital angular momenta then become constant. Therefore, the final 
semi-major axis for the HJ is w 2r p , where r p is the pericenter distance at the beginning of the 
capture phased. 

The same type of evolution shown in Figure 2 is seen with a broad range of initial conditions. 
There are two main routes to forming a HJ through the dynamical evolution of the systems we 
consider here. In the first, tidal friction slowly damps the growing eccentricity of the inner planet, 
resulting in circularized, prograde HJs. In the second, a sudden high-eccentricity spike in the 
orbital evolution of the inner planet is accompanied by a flip of its orbit. The planet is then quickly 
circularized into a retrograde short-period orbit. We can estimate the relative frequencies of these 
two types of outcomes using Monte Carlo simulations. Given the vast parameter space for initial 
conditions, a complete study of the statistics is beyond the scope of this Letter (but see Naoz et 
al., in preparation). However, we can provide a representative example: consider systems where 
the inner planet was formed in situ at ax — 5 AU with zero obliquity (orbit in the stellar equatorial 
plane) and with some small eccentricity e\ = 0.01, while the outer planet has 02 = 51 AU. The 
masses are m\ = 1 Mj and m 2 = 3 Mj. We draw the eccentricity of the outer orbit from a uniform 
distribution and the mutual inclination from a distribution uniform in cosi between and 1 (i.e., 
isotropic among prograde orbits). For this case we find that, among all HJs that are formed, about 
7% are in truly retrograde motion (i.e., with respect to the total angular momentum) and about 50% 
are orbiting counter to the stellar spin direction. Note that the latter fraction is significantly larger 
than what previous studies have obtained with stellar-mass perturbers (at most ~ 10%™). The 
high observed incidence of planets orbiting counter to the stellar spin direction 1 may suggest that 
planet-planet secular interactions are an important part of their dynamical history. 

Our mechanism requires that two coupled orbits start with a relatively high mutual inclination 
(i > 50°). The particular configuration in Figure 2 has a very wide outer orbit similar to those of 
directly imaged planets such as Fomalhaut bP^and HR 8799b 27 . In this case the inner Jupiter could 
have formed in its original location in accordance with the standard core accretion model^ on 
a nearly circular orbit. An alternative path to such a configuration involves strong planet-planet 
scattering in a closely packed initial system of several giant planets 7 . Independent of any particular 
planet formation mechanism, we predict that systems with misaligned HJs should also contain a 
much more distant massive planet or brown dwarf on an inclined orbit. 
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Figure 1: Dynamical evolution of a representative planet and brown dwarf system. Here we ignore 
tidal dissipation, but we do include the lowest-order post-Newtonian precession rate for the inner 
orbit. Here the star has mass 1 M & , the planet 1 Mj and the outer brown dwarf 40 Mj. The inner 
orbit has ai = 6 AU and the outer orbit has a 2 = 100 AU. The initial eccentricities are ei = 0.001 and 
e 2 = 0.6 and the initial relative inclination i = 65°. We show from top to bottom: (a) the inner orbit's 
inclination (h); (b) the eccentricity of the inner orbit (as 1 - ei); (c) and (d) the z-component of 
the inner- and outer-orbit's angular momentum, normalized to the total angular momentum (where 
the z-axis is defined to be along the total angular momentum). The thin horizontal line in (a) 
marks the 90° boundary, separating prograde and retrograde orbits. The initial mutual inclination 
of 65° corresponds to an inner and outer inclination with respect to the total angular momentum 
(parallel to z) of 64.7° and 0.3°, respectively. During the evolution, the eccentricity and inclination 
of the inner orbit oscillate, but, in contrast to what would be predicted from evolution equations 
truncated to quadrupole order [shown by the thin curves in panels (a) and (b)], the eccentricity of 
the inner orbit can occasionally reach extremely high values and its inclination can become higher 
then 90°. The outer orbit's inclination always remains near its initial value. We note that more 
compact systems usually do not exhibit the same kind of regular oscillations between retrograde 
and prograde orbits illustrated here, as chaotic effects become more important and are revealed 
at octupole order (see Fig. 2). We find that ~ 50% of the time the inner orbit is retrograde. 
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Figure 2: Dynamical evolution of a representative two-planet system with tidal dissipation included. 
The inner planet becomes retrograde at 1 12 Myr, and remains retrograde after circularizing into 
a HJ. Here the star has mass 1 M , the inner planet lMj and the outer planet 3Mj. The inner 
orbit has a\ = 6 All and the outer orbit has a 2 = 61 AU. The initial eccentricities are e\ = 0.01 and 
e 2 = 0.6, the initial relative inclination i = 71.5°, and argument of periapsis is 45°. We show: (a) 
the inner orbit's inclination (ix); (b) the eccentricity of the inner orbit (as 1 - ei); (c) the semi-major 
axis, peri-, and apo-center distances for the inner orbit and the peri- and apo-center distances for 
the outer orbit; (d) the magnitude of the angular momentum of the inner orbit; and, in (e) and (f) 
the z-components of the inner and outer orbit's angular momenta, normalized to the total angular 
momentum. The initial mutual inclination of 71.5° corresponds to inner- and outer-orbit inclinations 
of 64.7° and 6.8°, respectively. During each excursion to very high eccentricity for the inner orbit 
[marked with vertical lines in panels (b) and (c)], tidal dissipation becomes significant. Eventually 
the inner planet is tidally captured by the star and its orbit becomes decoupled from the outer 
body. After this point the orbital angular momenta remain nearly constant. The final semi-major 
axis for the inner planet is 0.022 AU, typical for a HJ. The thin curves in panels (a),(b),(d),(e) and 
(f) show the evolution in the quadrupole approximation (but including tidal friction), demonstrating 
that the octupole-order effects lead to a qualitatively different behavior. For the tidal evolution in 
this example we assume tidal quality factors Q± = 5.5 x 10 6 for the star and Qj = 5.8 x 10 6 for 
the HJ (see Supplementary Information). We monitor the pericenter distance of the inner planet 
to ensure that it always remains outside the Roche limit 29 . Here, as in Figure 1, we also include 
the lowest-order post-Newtonian precession rate for the inner orbit. 
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Supplementary Information 



Octupole-order Evolution Equations and Angular Momentum Conservation 

Our derivation corrects an error in previous Hamiltonian derivations of the secular evolution 
equations. 

We consider a hierarchical triple system consisting of an inner binary (mi and m 2 ) and a 
third body (m 3 ) in a wider exterior orbit. We describe the system using canonical variables, known 
as Delaunay's elements, which provide a particularly convenient dynamical description of our 
three-body system. The coordinates are chosen to be the mean anomalies, l\ and / 2 , the longitudes 
of ascending nodes, hi and h 2 , and the arguments of periastron, gi and g 2 , where subscripts 1, 2 
denote the inner and outer orbits, respectively. Their conjugate momenta are: 



Li = mim2 J k 2 {mi + m 2 )ai , (1) 
mi + m 2 



m 3 (mi + m 2 ) r— ■ ; — 

Jk 2 (m 1 +m 2 + m 3 )a 2 

mi + + ^3 



Gi = L^l-el, G 2 = L 2 ^l-e 2 2 , (2) 
where k 2 is the gravitational constant, and 

Hi — G\ cos i\ , H 2 = G 2 cosi 2 , (3) 

where Gi and G 2 are the absolute values of the angular momentum vectors (Gi and G2), and Hi 
and H 2 are the z-components of these vectors. 

We choose to work in a coordinate system where the total initial angular momentum of 
the system lies along the z axis. The transformation to this coordinate system is known as the 
elimination of the nodes 3017 ; the x-y plane in this coordinate system is known as the invariable 
plane. Figure[3]shows the resulting configuration of the orbits. We obtain simple relations between 
Hi, H 2 , d and G 2 , using Gtot — Gi + G 2 ~. 

C0S! = 2Gk ■ <4> 

/-i2 1 /~i2 /~i2 

H i = 7^ > (5) 

ri2 1 r*2 r*2 

H 2 = ^ , (O) 

^Wot 
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where the relation for Hi comes from setting G 2 = G tot — Gi (and similarly for H 2 ). Because total 
angular momentum is conserved by the evolution of the system, we must have Gi(t) + G 2 (t) = 
G to t = G tot z, implying that 

h 1 (t) = h 2 (t)-n. (7) 

\3tot 




The invariable plane 



Figure 3: The angular momenta of the bodies after the elimination of the nodes (see also Ref. 4). Note 
that all three vectors are in the same plane. The mutual inclination i = %i + % 2 is the angle between 
Gi and G 2 . 



The Hamiltonian for the three-body system can be transformed into the form 

H = U?(L 1 ) + H 2 X (L 2 )+H 12 , (8) 

where and K 2 represent the Keplerian interaction between bodies 1 and 2 and the central 
body, and T-L 12 represents the interaction between body 1 and body 2. The Kepler Hamiltonians 
depend only on the momenta L x and L 2 , while the interaction Hamiltonian, Hi 2 , depends on all 
the coordinates and momenta. Due to the rotational symmetry of the problem, %i 2 depends on hi 
and h 2 only through the combination hi — h 2 . Because we are interested in secular effects, we 
average the Hamiltonian over the coordinates (angles) li and l 2 , obtaining the secular Hamiltonian 

H = H?(L 1 )+H?(L 2 )+Hi2, (9) 

where 

#12 = ^ fj dh fj dl 2 Hi 2 . (10) 

For simplicity we first focus on the quadrupole approximation, where the error is more easily 
shown; it is then straightforward to see its effects at all orders in the hierarchical triple system's 
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secular dynamics expansion. The quadrupole Hamiltonian results from expanding H\ 2 to second 
ordei0in ai/a 2 : 

= + o(±)\ (ii) 

\a 2 J 

The resulting quadrupole-order Hamiltonian, ), depends only on the coordinates gi, hi, and h 2 , 
with the latter two appearing only in the combination hi — h 2 : 

n i S = nS{gx,h 1 -h 2 ). (12) 

Previous calculations 17 ' 20 eliminated hi and h 2 from the Hamiltonian using eq. ©, obtaining a 
quadrupole Hamiltonian that depends only on gi. But, this is incorrect! Such a Hamiltonian 
would imply that all quantities in eq. © are constant except G\, i.e. that eq. © is incorrect. Thus 
the previously used formalism did not conserve angular momentum. The initial Hamiltonian is 
spherically symmetric, and therefore does conserve angular momentum; the correct quadrupole 
Hamiltonian does as well. Because the correct quadrupole Hamiltonian depends on hi and h 2 
through the combination hi — h 2 , we have 

Hi = ~H 2 , (13) 

or 

Hi + H 2 = Gtot = const. (14) 



The mathematical error affects all orders in secular perturbations. The independence of the 
secular quadrupole Hamiltonian on h X2 was the source 17 of the famous relation cos 2*1,2 yl — ef 2 = 
const. In the correct derivation, this relation does not always hold. However, in a certain limit, it 
does. From eq. ©, we see that 

Hi = ^P-Gi - ^G 2 . (15) 
'j'tot ^tot 

When G 2 ~ G tot 3> G\, we have 

G 2 ■ 

Hi « --^G 2 . (16) 

At the quadrupole level is independent of g 2 , so G 2 = 0, implying 

Hi » 0, (17) 

when G 2 ~ G tot ^> Gi. This is precisely the limit considered in previous works 2,3 ' 17 ' 18 ' 19 ' 31 , so 
their conclusion that Hi^ 2 = cosii j2 Jl — ef 2 = const is correct (though not for the reason they 
claim), but the limit where G 2 ^> Gi is not sufficient for our work. 

In some later studies, the assumption that Hi = const was built into the calculations of 
secular evolution for various astrophysical systems 9,22-24 , even when the condition G 2 3> G\ was 

^The first order term in ai /a2 averages to zero, so the quadrupole term is the first term to contribute to H12 
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not satisfied. Moreover many previous studies simply set i 2 = 0, which is repeating the same 
error. In fact, given the mutual inclination i, the inner and outer inclinations i\ and i 2 are set by the 
conservation of total angular momentum: 



C0S2 2 = 



C0S2i = 




(18) 



(19) 



Tidal Friction 

We adopt the tidal evolution equations of Ref. 16, which are based on the equilibrium tide 
model of Ref. 32. The complete equations can be found in Ref. 2, eqs A1-A5. Following their 
approach (see their eq. A10) we set the tidal quality factors Q lj2 oc P in [see also Ref. 33]. This 
means that the viscous times of the star and planet remain constant; the representative values we 
adopt here are 5 yr for the star and 1.5 yr for the planet, which correspond to Q± = 5.5 x 10 6 and 
Qj = 5.8 x 10 6 , respectively, for a 1-day period. 

Comparison to Observations 

The observable parameter from the Rossiter-McLaughlin effect is the projected angle be- 
tween the star's spin and the orbital angular momentum (the projected obliquity) 14 . Here instead 
we focus on the true angle between the orbital angular momentum of the inner planet and the total 
angular momentum. Projection effects can cause these two quantities to differ in magnitude, or 
even sign. 

Moreover, several mechanisms have been proposed in the literature that could, under certain 
assumptions, directly affect the spin axis of the star. These mechanisms can re-align the stellar 
spin axis through tidal interactions with either a slowly spinning star 29 or with the outer convective 
layer of a sufficiently cold star 10 . Additionally, a magnetic interaction between the star and the 
protoplanetary disk could also lead to misalignment between the stellar spin and the disk 6 . 

These effects can potentially complicate the interpretation of any specific observation. Nev- 
ertheless, if hot Jupiters are produced by the simple mechanism described here, many of their orbits 
should indeed be observed with large projected obliquities. 
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